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Disks and Winds in Young Stars 


a) Wind-Disk Interactions 

Our paper on physical models for the winds of T Taun stars has been accepted by the 
Astrophysical Journal (Hartmann and Raymond 1989). 

b) Disks in T Tauri stars 

The low-mass, pre-main sequence T Tauri stars exhibit large excesses of infrared 
emission from warm dust grains (T — 100 K to 1500 K). This material must extend close to 
the star to explain the high dust temperatures and have a large radial extent to produce a 
range of temperatures, yet there must be relatively clear lines of sight to the stellar 
photospheres of these objects. A simple way to deal with the observations is to assume that 
the dust is in a disk around the young star rather than a series of spherical shells. Several 
other lines of evidence indicate that most young stars are surrounded by dusty disks, including 
direct observations of extended emission (Beckwith et al. 1984; Grasdalen et al. 1984). The 
general blue-shift of forbidden line emission in many young stars has been attributed to 
occultation of receding material by an opaque disk (Appenzeller, Jankovics, and Ostreicher 
1984), and bipolar flow models satisfactorily account for the observed line profiles (Edwards 
et al. 1987). 

Adams and Shu (1986) pointed out that a flat, dusty disk will absorb and reprocess - 
25% of the photons emitted by the central star, and Adams, Lada, and Shu (1987) later 
applied this idea to explain the infrared (IR) energy distributions of pre-main sequence stars. 
Adams, Lada, and Shu noted that the observed IR excesses are considerably larger than 25% 
of the stellar flux, suggesting that some of the excess might be produced by a flow of material 
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through the disk (i.e., accretion). 

We have analyzed spectral energy distributions of T Tauri stars to place limits on the 
amount of accretion which might occur during this early phase of stellar evolution (Kenyon 
and Hartmann 1987, KH). The typical system has a far-IR energy distribution which follows 
XF^ «»« X~° 7 , while the flat disk model predicts XF^ « X" 1-3 . If the height of the disk 
photosphere increases faster than linearly with radius (i.e., the disk “flares”), then the energy 
distribution of a typical T Tauri star will be flatter than XF^ X because disk annuli at 
large radii from the central star can absorb more photons. A physically plausible amount of 
flaring results in an energy distribution which matches the observations of the typical system. 

Possibly important uncertainties in the these calculations are the assumptions that the 
disk absorbs all of the photons incident upon it and radiates as a blackbody at the local 
equilibrium temperature. To address these issues, we have determined the temperature 
distribution of a dust layer in radiative equilibrium heated from above by radiation obliquely 
falling on the surface. We illuminated a plane-parallel slab of dust with radiation from a 
4000 K star at inclinations and distances appropriate for a flat dusty disk having an outer 
radius of 1000 times the stellar radius, R*. The disk was divided into 16 annuli equally 
spaced in log radius. 

The dust opacity used to calculate the emergent radiation field is an important ingredient 
for the radiative transfer models. The opacity depends on the composition and size 
distribution of grains, which are very uncertain in the conditions expected for a circumstellar 
disk. For our first calculation, we have assumed that the grains are similar in composition 
and optical properties to grains found in the interstellar medium. We adopted results 
summarized by Savage and Mathis (1979) for wavelengths below 9000 A and those of 
Hildebrand (1984) for X > 27 pm. At intermediate wavelengths, X - 1-27 pm, we used the 
results of Jones and Merrill (1976), scaled to Savage and Mathis at 9000 A. In this 
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calculation, the opacity does not depend on the disk temperature. 

Our initial results are summarized in Figure 1. The bottom panel of Figure 1 shows the 
radial temperature distribution expected for a blackbody disk (solid line; from KH) and that 
determined from the radiative transfer calculation (dotted line). The radiative transfer result is 
defined as T(x x = 2/3), where X is the wavelength at which F x is a maximum. It is apparent 
that the blackbody temperature is a good approximation for most radii, but begins to fail for 
log R/R* > 2. At these large disk radii, the expected temperature is such that the Planck 
function peaks at a wavelength where the disk is nearly optically thin. The disk does not emit 
much radiation at these long wavelengths, X > 60 |im, so the peak of the energy distribution 
shifts to shorter wavelengths. At these wavelengths, the disk is optically thick and the 
radiation temperature is higher then predicted by the simple models. 

The flux distribution of the star plus model disk is compared to the blackbody model in 
the top panel of Figure 1. The blackbody is a fairly good approximation, but it tends to 
produce more long wavelength flux than the radiative equilibrium calculation. The small 
excess of radiation relative to the blackbody calculation at 1 pm is a result of stellar flux 
scattered off the disk. At longer wavelengths, the new calculation has a deficiency of flux 
compared to the blackbody, because of optical depth effects. 

Our initial results suggest that energy distributions for circumstellar disks are reasonably 
well-approximated by the optically thick blackbodies used in our original calculations. In the 
next year, we plan to see how sensitive the results are to inclination angle in the disk and 
changes in the opacity law. These radiative transfer solutions will give us a better 
quantitative understanding of reprocessing disks surrounding T Tauri stars. Constraints on the 
amount of far-IR flux that can be produced by reprocessing disk models will provide 
important limits on the necessity for normal viscous accretion in the majority or TTS and for 
the non-viscous transport of mass and angular momentum as a solution to the peculiar far-IR 
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energy distributions of several T Tauri stars. 
c) T Tauri Stars in the HR Diagram 

The Hertzsprung-Russell diagram has been a basic tool for understanding of how single 
stars evolve once they are on the main sequence. Results for ITS are more ambiguous, 
because it has been difficult to estimate the luminosity and effective temperatures of the 
stellar photospheres in these objects. Most TTS display significant infrared and ultraviolet 
excess radiation over the energy produced by a normal stellar photosphere, and recent results 
suggest that this radiation is produced by a circumstellar disk which surrounds the star. This 
extra radiation source complicates estimates of the intrinsic stellar luminosity and the 
interpretation of the HR diagram. 

We have nearly completed a study to assess the uncertainties disks introduce in assigning 
luminosities and effective temperatures to TTS. An obvious problem with placing a disk 
around a star is that the disk occults a portion of the stellar hemisphere when the system is 
viewed at an inclination other than pole-on. Even if the disk did not emit any radiation itself, 
the apparent stellar luminosity, L app = 4 n d 2 F - where d is the distance and F is the 
bolometric flux received at earth, will be smaller than the true luminosity, L*, when the disk 
inclination is i * 0°. The addition of possible accretion luminosity complicates estimating the 
predicted distribution of L app assuming all TTS have the same L*. 

Our approach to this problem is to predict the distribution of L app assuming that all TTS 
have the same intrinsic luminosity, L»; TTS are randomly distributed in inclination, j; and 
TTS have some accretion luminosity, L^. We have constructed a sequence of model energy 
distributions for TTS, following methods outlined by KH. These models consider three 
radiation sources: a star having mass, M*, radius, R*, and temperature, T*; a steady-state disk, 
in which the accretion rate, M, is constant with radius through the disk; and a geometrically 
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thin boundary layer , where disk material rotating at Keplerian velocities loses kinetic energy 
and comes to rest on the slowly-rotating central star. 

Preliminary two-dimensional hydrodynamic calculations indicate that it is reasonable to 
consider the boundary layer as a ring of width /R* and scale height H bl (f « 1; H w « R*) 
with/= H bl /R*. We assume that the region can be characterized by a single effective 
temperature, T bl , and density, n bl . The boundary layer density formally is n bl - Z max /2H bl , but 

is sensitive to the choice of the viscosity parameter, a, and the mass accretion rate, M. 
We have chosen to determine the boundary layer structure from the energy conservation and 
continuity equations: 

-2^1 = 4 7c/R* J JtB x (T b i) [1 " e~ Tx * Thl) ] dX , (1) 

2R* 0 

n bl = M - - . (2) 

47t Vr p R* H b j 

where v R is the radial drift velocity through the boundary layer and 1*. is the optical depth. 

KH adopted v R = 1 km s _1 as a compromise between the local sound speed - 10 km s -1 and 
the drift velocity in the disk - 0.1 km s -1 , and this seems to be a reasonable assumption until 
the dependence of and other physical variables in the boundary layer on assumptions 
concerning the viscosity are understood. 

The “structure equations” for the boundary layer cannot be solved exactly for H bb n b j, 
and T w unless the region is optically thick at all wavelengths. We adopted initial guesses for 
these variables assuming the optically thick limit and applied a Newton-Raphson technique to 
determine self-consistent values. The procedure typically converges in 2-3 iterations and 
solves equation (1) to ~ 0.001%. 

As in KH, we assume that the disk can be considered as a collection of concentric, 
optically thick annuli, each of which radiates as a blackbody having an effective temperature, 
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T d . The flux emitted by each annulus is the sum of stellar energy absorbed and reprocessed 
by the annulus, Fr. and the energy liberated by viscous stresses within the disk (i.e., the 
accretion energy), F A . The effective temperature is computed according to the blackbody 
relation, T d 4 (R) = [F A (R) + F D (R)]/a. 

A sample model optical spectrum is compared to dereddened optical spectrophotometry 
of GG Tau in Figure 2. The model in the top panel of the Figure is for an accretion rate of 
M = 5.2 x 10 -8 Mq yr -1 onto a pre-main sequence star with a radius of 3 R©, a mass of 0.5 
Mo, and an effective temperature of 3750 K. The model has been scaled to a distance of 160 
pc to compare with the data for GG Tau in the bottom panel. The observations of GG Tau 
have been corrected for a modest amount of interstellar reddening (A v ~ 0.5 mag). 

It is apparent that the model adequately accounts for several features observed in GG 
Tau, including the height of the B aimer jump and the depths of several absorption bands. We 
have not tried to determine a “best-fit” to the data in the least-squares sense, but plan to 
develop software for this task over the next several months. 

To construct a theoretical H-R diagram for TTS, we have calculated spectral energy 
distributions for a variety of accretion rates M*M through a disk surrounding an MO pre-main 
sequence star with a radius of 3 R© and a luminosity of 1.6 L©. We characterize each 
calculation by the amount of “veiling” accretion produces at 5550 X, r, as a function of the 
inclination, i . The veiling parameter, r, is defined as the fraction of flux at 5550 X supplied 
by the disk and the boundary layer and ranges from 0 (no veiling) to 1 (complete veiling). 

The “stellar luminosity” of each model is estimated by normalizing the energy distribution of 
an M0 star to the calculated energy distribution at V, R, I, J, and K. This allows us to 
calculate the ratio of the inferred luminosity, L, to the true stellar luminosity, L*, at each 
wavelength, for a range of inclinations, /, and veiling, r. 
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If we assume that TTS are oriented randomly in space, then the probability density of 
TTS having log (L/L*) is: 

0.75 

p(log L/L*) = J n(r) A(r) — dr , (3) 

o 

where A(r) is a normalization constant. 

To compare the expression for p(log L/L*) with observations of TTS requires some 
knowledge of the errors introduced in estimating log (L/L*) from the data. Obvious sources 
of error include the distance, the interstellar extinction, and the spectral type, and our 
estimates for these are listed in Table 1. These results have been added in quadrature to 
produce a total error, o(log L/L*), which is listed in the final column of Table 1. The 
uncertainty in the luminosity is smallest in I and J, -0.1 in log (L/L*), and is severe for V, — 
0.2 in log (L/L*). 

Table 1 

Luminosity Dispersions for T Tauri Stars 


Wavelength 

Distance 

Extinction 

Spectral Type 

Total 

V 

0.058 

0.186 

0.064 

0.205 

I 

0.058 

0.090 

0.045 

0.116 

J 

0.058 

0.052 

0.073 

0.107 


We assume that the error in log (L/L*) is distributed randomly, and thus convolve the 
theoretical distribution for p(log L/L*) with a gaussian of width o(log L/L*) for comparison 
with the data. The final distribution is then: 

P(logL/L.) = J p(logl/l.) e H<10 ^- “ k>s ‘ JUy ° 1 ’) dlogl/1. , (4) 

Several examples of predicted probability densities are presented in Figure 3. The two 
lefthand panels show P(log L/L*) for models in which the veiling flux is 0 % (r = 0) and 50 % 
(r = 0.5) of the stellar luminosity at 5550 A. When the accretion luminosity is small (i.e.. 
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when r is close to zero), most of the disk luminosity is produced by reprocessing and the 
apparent stellar luminosity is close to the true stellar luminosity. The apparent luminosity 
greatly exceeds L* when accretion is a significant luminosity source, as shown in the upper 

left panel of Figure 3. 

It is not possible to predict the distribution of accretion rates, M, among TTS, so we 
have estimated the amount of veiling directly from observations of TTS. The observations are 
consistent with TTS having a roughly uniform distribution of veiling from r - 0 to r - 0.5 
with a few extra TTS having r - 0.5 to 0.75. We have considered two cases: (i) n(r) = 
constant from r = 0 to r = 0.5, and (ii) n(r) = 0.75 - r, for 0 < r < 0.75. Our results are 
displayed in the right panels of Figure 3, and it is obvious that the predicted probability 
distributions are not sensitive to the form of the veiling distribution. 

The models suggest that the combined effects of accretion and disk inclination would 
produce a measurable spread in the observed luminosities of TTS even if the underlying stars 
all had the same luminosity, L*. The width of the distribution is sensitive to the wavelength 
at which L* is estimated. We found that the minimum width occurs for the I and J 
bandpasses, so these wavelengths provide the best estimate for the stellar luminosity of a 
given TTS. Large dispersions in L/L* result at V and K, and these bandpasses should not be 
used to estimate the stellar luminosity of a TTS. 

The width of the probability distribution provides a useful estimate of the minimum age 
spread which can be measured in a T association. The luminosity of a star contracting to the 
main sequence along the Hayashi track follows L — 1 (t/3 x 10^ yr) ^ L© , so the minimum 
age spread which can be inferred from observations of L is ~ 1-2 x 10 6 yr. 

We have estimated the stellar luminosities, L(J), of TTS using published J magnitudes 
and extinctions, and have plotted the results as a histogram in Figure 4. The width of the 
observed distribution is significantly larger than that of the predicted distribution. An age 
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spread of - 3 x 10 6 yr is sufficient to explain the width of the observed distribution in L(J). 

We plan to complete this project over the next several months. These results have an 
obvious impact on the stellar birthline proposed by Stahler (1983) and one of our goals is to 
explore the statistical significance of the birthline for TTS. In addition, Walter et al. (1987) 
have discovered a class of pre-main sequence stars, the weak-lined TTS, which appear to lack 
(or possess minimal) ultraviolet and infrared excesses. Our analysis may make it possible to 
compare the luminosity functions for the classical and weak-lined TTS and determine if the 
weak-lined objects are TTS in which the disk has dispersed or whether some process has 
prevented the formation of disks around these stars. 
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Figure Captions 

Figure 1 - Comparison of predicted temperature (lower panel) and energy distributions (upper 
panel) for a blackbody disk (solid lines) and the radiative transfer calculation (dotted lines). 
The temperature at x = 1 in the transfer calculation follows the blackbody model until the disk 
becomes optically thin at far-IR wavelengths. 

Figure 2 - top panel: model energy distribution for a TTS with an optically thin boundary 
layer, bottom panel: observed energy distribution of GG Tau, a moderately veiled TTS. 

Figure 3 - Predicted luminosity functions for T Tauri stars. 

Figure 4 - Comparison of the observed distribution of TTS luminosities with predictions of 
the disk models describe din the text. 
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